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1  Introduction 


The  least  squares  method  plays  an  important  role  in  drawing  the  inferences  about  the  pa¬ 
rameters  in  the  nonlinear  regression  model.  In  this  paper  we  consider  the  least  squares 
estimators  (LSE’s)  of  the  following  sinusoidal  time  series  regression  model: 

Y {t)  =  Aocos{u>ot)  +  Bosin{uJot)  +  X {t}]  t  =  (1) 

Here  Aq  and  Bq  are  unknown  fixed  constants,  uq  is  an  unknown  frequency  lying  between  0 
and  TT.  X(t)’s  are  stationary  time  series  satisfying  the  following  assumption: 

Assumption  1 

£  oi{j)e{t-j),  £  |a0)|<oo  (2) 

j=z  — oo  j=-00 

where  e(t)’s  are  independent  and  identically  distributed  (i.i.d.)  random  variables  with  mean 
zero  and  finite  variance  >  0.  Here  ‘=’  means  X(t)  has  that  almost  sure  representation. 

This  is  an  important  and  well  studied  model  in  Time  Series  and  Signal  Processing  litera¬ 
ture.  See  for  example  Stoica  (1993)  for  an  extensive  list  of  references  for  different  estimation 
procudures.  Hannan  (1971,  1973),  Walker  (1969,  1971),  Kundu  (1993,  1995),  Kundu  and 
Mitra  (1995, 1996)  also  considered  this  or  similar  kind  of  model  to  study  the  asymptotic  prop¬ 
erties  of  the  different  estimators  and  some  of  the  computational  issues  have  been  discussed  in 
Rice  and  Rosenblatt  (1988).  Walker  (1971)  considered  the  approximate  least  squares  estima¬ 
tors  (ALSE’s)  and  proved  the  strong  consistency  and  the  asymptotic  normaliy  of  the  ALSE’s 
under  the  assumptions  that  the  errors  are  i.i.d.  random  variables  with  mean  zero  and  finite 
variance.  The  result  has  been  extended  by  Hannan  (1971,  1973)  to  the  case  when  the  errors 
are  stationary  random  variables  with  continuous  spectrum.  Kundu  (1993)  also  considered  a 
similar  model  and  proved  directly  the  consistency  and  the  asymptotic  normality  of  the  LSE’s 
under  the  assumption  that  X(t)’s  are  i.i.d.  with  mean  zero  and  finite  variance  and  they  are 
normally  distributed.  The  result  was  extended  to  the  case  of  general  mean  zero  and  finite 
variance  i.i.d.  errors  in  Kundu  and  Mitra  (1996).  In  this  paper  we  generalize  the  result  of 
Kundu  and  Mitra  (1996)  to  the  case  when  the  errors  are  coming  from  a  mean  zero  and  finite 
variance  stationary  process.  We  prove  directly  the  consistency  and  the  asymptotic  normality 
of  the  LSE’s  when  the  X(t)’s  satisfy  Assumption  1.  It  is  important  to  observe  that  we  do 
not  need  the  continuity  assumption  of  the  spectrum.  Our  approach  is  straight  forward  and 
different  from  that  of  Walker  (1969,  1971)  or  Hannan  (1971,  1973).  Hannan  (1971,  1973) 
obtained  the  result  for  the  one  parameter  case  after  making  the  Fourier  transform  of  the 
data.  We  observe  that  it  is  not  necessary  to  make  the  Fourier  transform  of  the  data.  We 
also  consider  the  multiparameter  case  and  obtained  the  explicit  expression  of  the  asymptotic 
covariance  matrix,  which  is  not  available  in  the  literature.  We  also  perform  some  numerical 
experiments  to  compare  the  small  sample  behavior  of  the  ALSE’s  and  the  exact  LSE’s.  In 
this  paper  the  almost  sure  convergence  means  with  respect  to  the  usual  Lebesgue  measure 
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and  it  will  be  denoted  by  a.s..  Also  the  notation  a  =  0{N^)  means 
N. 


is  bounded  for  all 


7^1 


The  rest  of  the  paper  is  organised  as  follows,  in  Section  2  we  prove  the  consistency  of 
the  LSE’s  and  establish  the  asymptotic  normality  results  in  Section  3.  The  results  for  the 
several  Harmonic  case  are  obtained  in  Section  4.  Some  numerical  results  are  presented  in 
Section  5  and  finally  we  draw  conclusions  in  Section  6. 


2  Consistency  of  the  LSE’s 


Let’s  denote  9n  =  to  be  the  LSE  of  9o  =  {Aq,Bo,uo),  obtained  by  minimizing 

N 

Q[^{9)  =  —  Acos{ut)  —  Bsin{uit))^  (3) 

t=i 

with  respect  to  9  =  (A,  B,  u).  It  is  important  to  observe  that  the  existence  and  the  uniqueness 
of  a  respective  measurable  function  satisfying  (3)  follows  along  the  same  line  of  Jennrich 
(1969).  To  prove  the  consistency  results  we  need  the  following  lemma. 


Lemma  1:  Let  X(t)  be  a  stationary  sequence  which  satisfies  Assumption  1,  then 


lim  sup 

N-yoo  Q 


t=i 


=  0 


a.s. 


(4) 


Before  giving  the  proof  in  details,  we  would  like  to  give  a  sketch  of  the  main  idea.  First  we 
show  that  (4)  holds  for  the  subsequence  N^.  Then  we  show  that 


sup  sup 

0  N^<K<{N+l)3 


^'^X{t)cos{t9) 

t-l 


i|;x(i)cos(w) 

t=i 


(5) 


converges  to  zero  a.s.  as  N  tends  to  oo. 


Proof  of  Lemma  1: 

-j  TV  1  OO 

—  f^X{t)cos{t9)  =  a{j)e{t  -  j)cos{t9)  = 

i^\  j=~oo 

Q:(7Xt  -  i){cos((t  -  j)9)cos{j9)  -  sin{{t  -  j)9)sin{j9)}  = 

Y  Oi{j)cos{j9)  Y  f  -  j)cos{{t  -  •j)9)  - 
4  Y  oc{j)sin{j9)Y^i^- j)9)  (6) 

■'''  j=-oo  t=l 
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Therefore 


1  oo  N 

Z  «O')cos(i0)X^e(f-i)cos((i-i)0)|  + 

0  -iV 


j  =—oo 
1  ^ 


tz=l 

N 


'iV  . 

J  =  — oo 


(7) 


e=i 


We  would  like  to  prove  that  both  the  terms  on  the  right  hand  side  of  (7)  converges  to  zero 
as  N  tends  to  infinity.  Now  observe  that 


1 


^sup|^  H  a{j)cos{j9)J2^{t-j)cos{{t-j)e)\ 


j=-oo 


N 

c 

t-l 


< 


-I  oo  f  N  'j  2 

^  XI  l«(i)n^sup|  j^e(t-i)cos((t- < 

^  j=-oo  I  t=l  ) 


iv.. 

1 


N 


E  |Q;(i)IS^+  E  E[\i:m€{Tn)e{Tn  +  t)\Y 

j=-oo  [  t=-N+l 


(8) 


where  the  sum  omits  the  term  t  =  0  and  the  term  is  over  A^— lt|  term  (dependent 

on  j).  Since 

X:  f;[|E^e(m)e(m  +  t)l]<  E  [|S^6(m)e(m  +  t)p]  ^  =  O(iVt)  (9) 

t=-N+l  t=-N+l 

(uniformly  in  j)  therefore  (8)  is  0{N~^).  Let  M  =  N^.  Therefore 


1 


N 


^sup|—  XI  a{j)cos{j9)Y,cos{{t-j)e)\'^  =  0{M  3) 

®  ^  j-~oo  f=l 

Similarly  the  result  is  true  if  the  cosine  function  is  replaced  by  the  sine  also.  Therefore 


(10) 


sup 

0 


1  M 

—  Y^X{t)cos{t9) 

t=l 


when  M  =  N^.  Now 


sup  sup 

0  N'^<K<{N+\)^ 


0  a.s. 


-E  ^  X{t)cos{te)  -  E  ^{'t)cos{t9) 


(11) 


K 


sup  sup  I E  ^ {i)cos{t9)  -  E  ^ {t)cos{t9)  + 

9  N^<K<{N+\)^  A/ 
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±f:Mt)cosm  -  ^■tx(t)cosm  < 

^  t=l 


E  W‘)l+  E  I^WI 


t=N3+l 


iV3  (AT  +  1): 


The  mean  squared  of  the  first  quantity  on  the  right  hand  side  of  (12)  is  dominated  by 
+  1)3  _  ;Y3j2  ^  0(Ar-2).  Similarly  the  mean  squared  of  the  second  quantity  on  the 

right  hand  side  of  (12)  is  dominated  by  =  0{N-^).  Therefore  both  will  converge  to 
zero  almost  surely,  which  proves  the  lemma. 

Corollary  1:  The  result  is  true  if  the  cosine  function  is  replaced  by  the  sine  function. 

Corollary  2:  It  can  be  proved  similarly  that  if  X(t)  is  a  sequence  which  satisfies  Assumption 
1,  then 

lim  sup  -^y]t^X{t)cos{t9)  =0  a.s.  (13) 

N-^OO  Q  JSI'^ 


Now  consider 


J^[QN{e)-QN{9o)]  = 

^  {(y(t)  -  Acos(a;f)  -  Bsin{ut))^  -  = 

^  t=i 

1  V  2 

—  Yi  iAocos{uJot)  +  Bosin{uJot)  -  Acos{ujt)  -  Bsin{ut))  + 

^  f=i 

2  ^ 

—  Y  (Aocos(u;ot)  +  BosiTi{u)ot)  —  Acos{ujt)  —  Bsin{u)t))  = 

^  t=i 

fN{A,  B,u))  +  (Jm{A,  B,uj) 


Now  with  the  help  of  lemma  1,  we  can  ecisily  conclude  that 

lim  sup  <7Ar(A,  BjCj)  =  0  a.s. 

where  the  set  Ss,m  for  (i  >  0,  is  as  follows; 


therefore  for  all  5  >  0 


lor  0  >  u,  IS  as  loiiows; 

SsM  =  {(^B,a;),lA-Ao|>5,|Al<M,|B|<M 

or  \B  -  Bo\>6,\A\<  M,\B\<  M 
or  [a;  —  cuqI  >  6,  |A1  <  M,  |B|  <  M} 


lim  inf  — 
Ss,M  N 


^  [Qn(^)  -  Q/v(0o)]  =  lim  sup  /;v(^,-B,a;)  >  0. 
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(17)  follows  easily  from  Kundu  and  Mitra  (1996).  Here  lim  means  limit  infimum.  Now 
suppose  be  the  LSE’s  of  {Ao,Bq,ujq)  and  they  are  not  consistent.  Therefore 

either 

Case  I;  For  all  subsequences  {Nk}  of  {N},  +  \Bnk\  tends  to  infinity  or 

Case  II:  There  exists  a  5  >  0  and  a  M  <  oo  and  a  subsequence  {Nk  }  such  that  {A^^,  B^^, 
^Nk)  S  ^5,Mi  tbr  all  K  1,  2,  . . .,. 

Now 

Qnk{-^Nk^^Nk^'^Nk)  ~  Qnk{^o,  Bo,uJo)  <0  (18) 

as  (^Nk)  Bis[^,u}nk)  is  the  LSE  of  (/Iq,  Bq,ujq),  when  N  =  Nk-  Observe  that  as  K  oo,  for 
both  the  cases,  the  left  hand  side  of  (18)  converges  to  a  number  which  is  strictly  positive, 
that  is  a  contradiction.  Therefore  the  LSE’s  of  the  model  (1)  have  to  be  strongly  consistent. 
Therefore  we  can  state  the  following  theorem: 

Theorem  1  If  9^  =  (Ai^,Bfj,Cjt>i)  is  the  LSE  of  the  nonlinear  regression  model  (1),  then 
it  is  a  strongly  consistent  estimator  of  6o  =  {Aq,  Bo,uo)- 


3  Asymptotic  Normality 


Q’nW 


(19) 


In  this  section  we  prove  the  asymptotic  normality  of  9i^  by  using  the  Taylor  series  expansion. 
Let’s  denote 

(SQn{9)  SQn{9)  6Qr,{e)\ 

[  6A  '  5B  '  6u  ) 
and  Q'li{9)  to  be  the  corresponding  3x3  matrix  which  contains  the  double  derivative  of 
Qn{9).  Therefore 

QV(»/»)  -  Q'M)  =  (e-  WM  (20) 

where  9  =  {A,  B,u)  is  a  point  in  the  line  joining  ^/v  and  9q.  Observe  that  although  9  depends 
on  N,  we  omit  it  for  brevity.  Since  Q'[^{9n)  =  0,  (20)  implies 


Now 


(»-e„)  =  -Q;(«„)[(3';(«) 


^%j^  =  -2f;.Y(()cosM) 

t-l 

SQn{9o) 

— —  =  -2  A  {t)sm{ujot) 
^Q/v(9o) 


t=i 

N 


5uj 


—  —2^tX{t){Aosin{uJot)  —  Bocos{ujot)) 


(21) 

(22) 

(23) 

(24) 


t=i 
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Also 


s^QN{e) 

5A^ 

6'^Qn{9) 

buP’ 


biobA 


S^QnW 

bubB 


b'^QN{0) 

bAbB 

Let’s  define 


2/  N  ^^Qn{0)  .a  .  2/  X 
=  25^  cos^(a;i),  =  2  snP{ujt), 

t=\  t=l 

N 

=  2  ^  t^[(AoCos(a;ot)  +  Bosin{uot)  —  Acos{u)t)  —  Bsin{u}t)  +  X{t))  x 

t=i 

{Acos{Cijt)  +  Bsin{u)t))  +  {Asin{ut)  —  Bcos{Cbt)Y]  (25) 

N 

=  2^t[sin{u)t){Aocos{u!ot)  +  BQsin{uot)  —  Acos{u)t)  —  Bsin{uit)  +  A’(t))  — 

t=i 

cos(Cjt){Asin{u}t)  —  Bcos{uit))\  (26) 

N 

=  —2  tlcos{Qt){Aocos{ujot)  +  Bosin{u>ot)  —  Acos{Qt)  —  Bsin{Qt)  +  A’(t))  — 

t=i 

sin{Cjt){Asin{u)t)  —  Bcos{Qt))]  (27) 

N 

=  2  ^  sin{u)t)cos{Gjt)  (28) 

t=i 


1  /V  j 

(Til  =  liin  —  cos^iuot)  =  - 

N^oo  N  ^  ^  2 

1  ^  1 

(T22  =  lim  t7  sirPiunt)  =  - 

JV-40O  yv  ^  ^  ’  2 

1  ^  1 

(T33  =  lim  TTT  X]  {Aosin{ujQt)  -  Bocos{LOot)f  =  -{AI  +  Bq) 

N-¥oc  “  b 

1  ^  1 
<7i3  =  (^31  =  ^  Y1  Botcos'^iujot)  =  -Bo 

1  ^  1 
^■23  =  0-32  =  -  lim  —  y]  Aotsin^iuot)  =  --Aq 

N^oo  4 

1  ^ 

(Ti2  =  (T21  =  lim  —  y^  sin(ujot)cos(u)ot)  =  0 
N-ioo  N  tP, 


Now  observe  that  as  a>  — )•  wq,  A  — >  Aq  and  B  Bq  a.s.,  we  have 


1  ^  _  1  ,  1 
lim  —  y^  coP(u)t)  =  lim  —  y^  cos'^ii^ot)  =  - 
/v-^oo  N  ^  ^  \r  ^  '  0 


N 


t=i  ^  t=i 

N  1  N 


1  _  1  1 
lim  —  y'  siTp(u)t)  =  lim  tt  y^  sin^iuot)  =  - 


(29) 

(30) 
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= Ji?.o -  Bcos{ot)y = 


1  ^  1 
Jin^  ^  S  {^osin{uJot)  -  BQCOs{ujQt)f  =  -(/ig  +  bI) 

/)!,“  W  %  [Asin(u,t)  -  Bcos(u,t))  = 

1  ^  1 

xl^”^  TH'il^BQCOS^{UQt)  =  -Bq 
N-yoo  4 

lim  :rT7T7— — r;::—  =  —  lim  -r;^y^ tsin(ujt)  (Asin(u)t)  —  Bcosiut))  = 
N-foo  2N^  5uj5B  n^oc  ^ 

1  ^  1 
-  E  tAosin^iuot)  =  --Ao 

1 


lim 


N-yoo  N  5A5B  N-yoo  N 
2  ^ 

lim  —  sm(a;oi)cos(a;oi)  =  0 
^-+oo  iV  “ 


2  ^ 

=  lim  —'^sin{iL>t)cos{uit)  = 


(31) 


(32) 


(33) 


(34) 


Let’s  define  the  3x3  matrix  S  =  {{cTij))',  i,j  =  1,2,3  and  also  define  the  3x3  diagonal  matrix 
D  as  follows  D  =  diag  {  }.  Rewrite  (21)  as 

{9  -  9o)D-^  =  -Q'jv(^o)D  [dQ';(9)D]~'  (35) 

Now  from  (29)  -  (34)  we  obtain 


lim  DQ'U9)D  =  lim  DQ'U9o)D  =  2S 

N-aoo  N-¥oo 


where 


Li 


-1^0  \{Al  +  Bl) 


and  E  ^  exists  if  (^Iq  +  -So)  >  0  il-  is  as  follows; 


E"^  = 


2^0  +  25o 

4 

1 

(OICaS 

o 

-3JB, 

—  ^AqBq 

IflJ  +  2Al 

3j4o 

-3Bo 

CO 

o 

6 

(36) 


(37) 


(38) 
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Now  from  the  Central  Limit  theorem  of  Stochastic  Process  (see  Fuller;  1976)  it  easily  follows 
that  (5)v(^o)D  tends  to  a  multivariate  (3- variate)  normal  distribution  eis  given  below; 

Q'M)^  ^  N^{0,4a^cE)  (39) 


where 

00  OO 

c  =  |  S  Q!(i)cos(a;oi)P  + 1  ^  a{j)sin{uoj)\'^ 

j=-oo  j^-OO 

Therefore  we  have; 

{§N  -  9o)D-^  ^  N^{0,a^c^-^)  (40) 

Now  we  can  state  the  result  as  the  following  theorem; 


Theorem  2  Under  the  assumptions  of  Theorem  1,  [N'i  —  ^o)  (J5/v  —  Bq),  {Cj^  — 

Wo)}  converges  in  distribution  to  a  3-variate  normal  distribution  with  mean  vector  zero  and 
the  dispersion  matrix  is  given  by  c  where  c  and  are  as  defined  before. 


4  Multiparameter  Case 


In  this  section  we  will  extend  the  results  of  Section  2  and  Section  3  to  the  following  model: 

M 

T(t)  =  y~^  Aficos(uJnt)  4-  Bn  sinju^i)  +  X(t);  t  =  (41) 

K=\ 

where  Aq,  Bq  are  arbitrary  real  numbers  and  Wq^’s  are  the  distinct  frqucncies  lying  between 
0  and  TT  for  K  =  1,  . . M.  X(t)’s  satisfy  Assumption  1. 

Let  us  use  the  following  notations  A  =  (A‘,...,A^),  B  =  (B^,...,B^)  and  u  = 
(w\...,w^).  Similarly  Aq,  Bq,  wq  and  An,  Bn  and  wn  are  also  defined.  We  would 
like  to  investigate  the  consistency  and  the  asymptotic  normality  properties  of  the  LSE’s 
obtained  by  minimizing  Rn{^)  =, 

N  /  M 

^  I  T(t)  —  ^  [A^cos{LO^t)  +  B^ sin{u>^^t)] 

t=l  V  K=l 

with  respect  to  $  =  (A,B,w).  Now  we  have  the  following  result: 


Theorem  3  =  (An,Bn,u)i\[)  is  the  LSE  of  =  (Ao^BqjUq),  then  $Ar  *5  a  strongly 

consistent  estimator  o/'I'o. 
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Proof  of  Theorem  3:  With  the  help  of  Lemma  1  and  using  the  similar  kind  of  techniques 
as  that  of  (Kundu  and  Mitra;  1995),  the  results  can  be  established. 


Let’s  denote  the  1  x  3M  vector  R^i^)  as  follows: 

(SR„(i)  SRy{'i)\ 

-sIT-  j 

and  R%{^)  denotes  the  3M  x  3M  matrix  which  contains  the  double  derivative  of  Rff{^). 
Now  we  have 

=  {^N  -  ^o)Rn{^)  (42) 

where  $  =  (A,  B,cu)  is  a  point  in  the  line  joining  and  $o-  Since  =  0,  we  have 

-  $o)  =  -R'^i^o)  (43) 

Let’s  define  the  3M  x  3M  diagonal  matrix  V  whose  first  2M  diagonal  elements  are  N~^  and 
the  last  M  diagonal  elements  are  N~'i.  Therefore  we  can  write  (43)  as 

($;v  -  $o)V-i  =  -i?;($o)V-'  [V-iii!';(l>)V-']"' 

Now  using  the  similar  kind  of  arguments  as  of  Section  3,  we  can  say  that 

ii!)v($o)V->  Asm  (0, 4(7^0 


where  G  is  a  3M  x  3M  matrix  and  it  has  the  following  structure 


Gii 

Gi2 

Gi3 

G21 

G22 

G23 

G31 

G32 

G33 

where  each  of  the  Gy  is  a  M  x  M  matrix  and 


Gi3  =  Gsi  =  diag 

G23  =  G32  = -dzar/ j-ToCi, . . . , 1 
G33  =  \diag{di,...,dM} 

D 

Gi2  =  0. 


(44) 


(45) 
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here  ck  = 


I  51  oc{j)cos{oj^ j)\‘^  +  \  5^  a(j)sm(w^i)P 

j=-oo  j=-oo 

and  dft-  =  for  K  =  1,  . . .,  M.  Observe  that 


lim  VR'LmV  =  lim  Vi2';(^>o)V  =  2r 

N^oo  '  N-^oo  ' 


here  the  3M  x  3M  matrix  F  is 


r  = 


jIm  0  Si 
0  jIm  S2 
Si  S2  S3 


where  Si,S2,S3  are  M  x  M  diagonal  matrices  as  follows; 

51  =  -^diag{Bi . Bo"} 

52  =  -jdiag{Ai,...,A^} 


S3  ^diug  {dj,  •  ■  • )  d^} 


and  Im  is  the  identity  matrix  of  order  M.  Since 

r~^  =  4 

where 


5R4  +  2Ri 

R2 

Ra 

R2 

iRi  +  2R4 

Rs 

Rs 

Rs 

Re 

—Zdiag 


Ri  =  diag 
R2 
R3 

R4 

Rs 
Re 


(BJ) 


1\2 


^diag 


diag  | 
3diag 


d,  ’■■■’  dM  j 

Ai^\ 

1  d,  ’•••’  dM  J 
lEi 

Idi”--’  dM  J 

di  ’■■■’  d 


M 


/do  <;] 

Idr’--’  dM  J 

6dmf/ [ 

Idi  dw  J 


we  have 


(cl.;v  -  <Po)V-'  ^  iV3M  (o,a2r-‘Gr-i) 
therefore  we  can  state  the  result  as  the  following  theorem; 


Theorem  4  Under  the  assumptions  of  Theorem  3,  {N^{Ai^  —  Ao)  ,A^5(Ba^  — Bo),  Ni{uN- 
(jJo)}  converges  in  distribution  to  a  3M-variate  normal  distribution  with  mean  vector  zero  and 
the  dispersion  matrix  is  given  by 


5  Numerical  Experiments 


In  this  section  we  perform  some  Monte  Carlo  simulations  to  see  how  the  asymptotic  results 
work  for  small  sample.  We  considered  the  following  model: 

Y{t)  =  Aocos{ujot)  +  Bosin{uot)  +  X{t)-,  t  =  l,...,N.  (51) 

We  took  Ao  =  Bo  =  1.5,  a;o  =  .25  tt  («  0.735398),  .50  tt  («  1.570796)  and  .75  tt  («  2.356194). 
X(t)  =  e(t)  +  .5e(t-l),  where  e(t)’s  are  i.i.d.  normal  random  variables  with  mean  zero  and 
variance  one.  Numerical  results  are  reported  for  N  =  10,  15,  25.  All  these  computations 
were  performed  at  the  Indian  Instituet  of  Technology  Kanpur,  using  PC-486  and  the  random 
deviate  generator  proposed  by  Press  et  al.  (1992).  For  a  particular  N  and  cj,  1000  different 
data  sets  were  generated  and  for  each  data  set  we  estimated  the  nonlinear  parameters  by 
two  different  methods,  one  (denoted  by  L.S.)  by  directly  minimizing  (3)  with  respect  to 
the  different  parameters  and  the  other  one  (denoted  by  A. L.S.)  by  first  making  the  Fourier 
transform  of  the  data  as  suggested  by  Hannan  (1971,  1973),  Walker  (1971).  We  computed  the 
average  estimates  and  the  average  mean  squared  errors  over  1000  replications.  We  reported 
the  result  in  Table  1  for  the  frequency  only  because  the  others  are  quite  similar  in  nature. 
The  figures  in  the  top  denote  the  average  estimates  and  the  figures  in  the  parenthesis  below 
give  the  corresponding  averge  mean  squared  errors.  We  also  computed  the  95%  confidence 
interval  for  uj  for  each  data  sets.  The  results  are  reported  in  Table  2.  The  first  figure  in 
the  parenthesis  is  the  average  length  of  the  confidence  interval  and  the  second  figure  is  the 
coverage  frequency  over  1000  replications.  From  Table  1  and  Table  2,  it  is  clear  that  although 
asymptotically  both  the  methods  are  same  but  for  small  sample  it  is  observed  that  the  exact 
LSE’s  are  better  than  the  ALSE’s.  The  average  mean  squared  errors  of  u  are  lower  for  the 
usual  LSE’s  for  almost  all  the  sample  sizes  and  for  all  u's.  About  the  confidence  intervals, 
it  is  observed  that  for  higher  values  of  ca,  the  confidence  interval  of  u  obtained  by  using  the 
exact  LSE’s  usually  give  higher  coverage  probability.  It  is  also  observed  that  for  both  the 
methods  as  N  increases  the  average  length  decreases  and  the  coverage  probability  increases. 


6  Conclusions 


In  this  paper  we  considered  the  one  parameter  and  multiparameter  sinusoidal  model  under 
the  assumption  of  additive  stationary  errors.  We  obtained  the  asymptotic  properties  of  the 
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LSE’s  directly  without  making  the  Fourier  transform  of  the  data.  We  also  obtained  the  ex¬ 
plicit  expression  of  the  covariance  matrix  for  the  multiparameter  case,  which  is  not  available 
in  the  literature.  From  the  numerical  study  it  is  observed  that  although  asymptotically  the 
two  mehtods  are  equivalent  but  the  exact  LSE’s  are  better  than  the  ALSE’s  in  terms  of 
the  mean  squared  errors.  Since  both  the  methods  require  the  same  amount  of  computaions, 
therefore  it  is  recommended  not  to  Fourier  transform  the  data  at  least  for  small  samples  to 
make  any  finite  sample  inference  from  the  asymptotic  results. 
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Table  1 


u 

N  =  10 

N  =  15 

N  =  25 

L.S.  A.L.S. 

L.S  A.L.S. 

L.S.  A.L.S. 

.257r 

.7093  .6949 

(.1314)  (.1367) 

.7525  .7190 
(.0287)  (.0539) 

.7871  .7794 

(.0078)  (.0139) 

.507r 

1.3402  1.2555 
(.3387)  (.5072) 

1.4371  1.4543 

(.1287)  (.1467) 

1.4749  1.4497 

(.0975)  (.1436) 

.757r 

1.7772  1.6292 

(1.060)  (1.426) 

2.1143  2.0450 

(.4455)  (.6202) 

2.2501  2.1875 

(.1987)  (.3637) 

Table  2 


UJ 

N  =  10 

N  =  15 

N  =  25 

L.S.  A.L.S. 

L.S  A.L.S. 

L.S.  A.L.S. 

.257r 

.507r 

.757r 

(.24,  46)  (.24,  .41) 
(.31,  53)  (.21,  .42) 
(.35,  62)  (.37,  .49) 

(.25,  .73),  (.15,  .53) 
(.29,  .78),  (.19,  .71) 
(.32,  .86),  (.21,  .80) 

(.15,  .88),  (.09,  .81) 
(.16,  .87),  (.10,  .73) 
(.17,  .94),  (.11,  .86) 
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